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ABSTRACT 


Previous  efforts  to  derive  body  motion  parameters  from  monostatic 
scattering  matrix  data  have  been  limited  by  mathematical  complexity. 

An  account  of  these  efforts  is  given.  Then  a  proposed  method  is 
described  for  the  case  of  precession,  which  allows  the  simplicity  of 
a  geometric  approach  and  avoids  approximation  in  its  concept. 

This  method  is  formulated  as  a  minimization  problem  with  only 
three  of  the  precession  parameters  as  the  arguments.  The  basic  idea 
is  to  find  a  cone  on  which  the  spacing  of  successive  positions  of 
the  body  axis  shall  be  the  dynamically  correct  spacing.  Then  the 
two  remaining  parameters  are  easily  evaluated. 
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SECTION  I 


INTRODUCTION 

In  the  course  of  seeking  methods  of  deriving  motion  parameters 
from  scattering  matrix  data,  a  number  of  methods  have  been  considered 
and  reported  earlier  by  the  authors  in  a  number  of  internal  MITRE  reports. 
These  methods  appear  to  be  limited  in  practice  either  by  computational 
complexity  or  by  the  need  for  approximation. 

In  view  of  such  difficulties,  it  is  expedient  to  consider,  at 
least  initially,  a  method  which  allows  the  simplicities  of  a  geometric 
approach  and  conceptually  avoids  approximation.  Since  it  is  a 
monostatic  method  without  approximation,  it  requires  the  assumption 
of  a  motion  model  which  in  this  case  is  precession  for  objects  of 
rotational  symmetry.  The  formulation  given  here,  however,  can  be 
made  in  terms  of  just  three  motion  precession  parameters.  In  this 
work,  the  assumption  is  made,  of  course,  that  such  an  object  depolarizes 
at  the  operating  frequency  to  an  extent  which  allows  the  determination 
of  the  orientation  of  the  body  (electrical)  axis  in  a  plane  normal 
to  the  radar  line  of  sight.  It  is  to  be  noted  also  that  a  monostatic 
method  for  determining  the  body  axis  orientation  using  any  type  of 
observation  data  for  a  body  of  axial  symmetry  yields  ambiguities. 

Before  discussing  the  method  itself,  it  is  useful  to  put  into 
perspective  the  general  problem  in  terms  of  the  motion  model. 
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observations,  unknowns,  and  analytic  formulation.  Thus,  the  paper 
divides  into  two  parts  with  background  and  alternatives  first  discussed 
followed  by  an  exposition  of  the  proposed  method. 

Finally,  it  is  to  be  noted  that,  although  we  deal  with  one 
motion  technique  applied  to  scattering  matrix  data,  the  use  of  phase 
and  polarization  information  in  radar  returns  has  found  widespread 
application  in  work  at  MITRE,  benefiting  methods  in  motion  and  shape 
determination  both  for  long  and  short  pulse  data  on  whole  object  or 

r  1  2I 

scattering  center  descriptions  ^  .  One  author  (0)  has  applied 

the  scattering  matrix  to  wedge  parameters  in  other  internal  MITRE 
reports . 
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SECTION  II 


THE  PROBLEM 


Using  an  earth  based  radar,  observations  are  made  of  an  object 
in  orbit  possibly  undergoing  additional  rotation  movement  about  its 
center  of  mass.  We  illustrate  the  three  motions  in  Figure  1. 

In  general  all  three  motions  act  to  determine  the  object’s 
orientation  with  respect  to  the  radar  so  that  this  target-sensor 
orientation  affects  the  various  possible  sensor  measurements.  In 
this  paper  we  deal  with  the  source  data  in  the  form  of  the  polarization 
scattering  matrix  from  a  monostatic  radar.  Written  in  a  circular 
polarization  basis,  this  matrix  can  be  represented  for  backscatter 
from  a  reciprocal  scatterer  as: 


(1) 


where  a,  3,  v  are  complex,  §  represents  body  orientation  about  the 


an  object  with  a  plane  of  symmetry  containing  the  line  of  sight  , 
a  =  3.  Use  of  monostatic  data  implies  both  a  n  and  TT/2  arioiguity 
in  §  for  axially  symmetric  bodies.  Since  the  sense  of  the  body 
axis  vector  is  not  of  interest  here,  the  TT  ambiguity  need  not  be 
dealt  with.  The  process  of  obtaining  a  solution  provides  a  means 
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in  the  monostatic  case  for  resolving  the  TT/2  ambiguity.  The  condition 
of  sufficient  depolarization  is  imposed  to  have  a  measurable.  This 
latter  condition  is  enhanced  at  operating  frequencies  near  resonance. 

As  the  frequency  increases  with  geometrical  optics  applying,  all 
targets  become  sphere-like  with  (X  0.  Finally  the  assumption  of 
backscatterer  is  also  an  approximation,  good  for  many  motion  environments. 

We  can  consider  the  motion  problem  using  scattering  matrix  data 
in  terms  of  two  planes  designated  the  viewing  plane  and  the  measuring 
plane  as  shown  in  Figure  2 . 

Considering  all  motion  effects  absorbed  into  the  measurement 
of  ? (t) ,  we  find 


<  ,  u(t)  > 

Y(t)  =  tan  ?(t)  =  -  (2) 

<  X^(t),  u(t)  > 


in  which  <,  >  means  scalar  product.  The  projections  are  indicated 

A  A 

in  Equation  (2)  to  be  onto  time  varying  vectors  y„(t)  and  ^  (t) . 

R  R 

This  time  dependence  implies  that  a  fixed-in-space  inertial  frame 
is  to  be  used  as  a  reference  system. 

Equation  (2)  may  be  recast  into  the  following  useful  forms 

<  CY(t)  >  =»  <  n(t)  ,  u(t)  >  =  0  (3) 
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Figure  2.  Motion  Using  Scattering  Matrix  Data 


in  which  n(t)  lies  in  the  measuring  plane  and  is  normal  to  the  viewing 

A 

plane.  With  u  taken  as  unrestricted  in  orientation  and  magnitude, 
Equation  (3)  is  the  equation  of  the  viewing  plane. 

Before  proceeding  to  more  detailed  analytical  formulations,  it 

is  worthwhile  to  view  the  motion  of  the  viewing  plane  as  decomposed 

A  A 

into  infinitesimal  rotations  of  u  and  r  which  are  unit  vectors  in 
the  direction  of  the  body  axis,  u,  and  line  of  sight,  r,  respectively. 
We  illustrate  the  situation  in  Figure  3. 

The  infinitesimal  rotation  of  the  viewing  plane  can  be  given 
in  terms  of  R  as 


dR  =  R  X  dH 


r  X  dn  +  u  X 


(4) 


*  dr  +  du 


Since  any  component  of  dO  normal  to  the  viewing  plane  keeps 

AAA 

this  plane  invariant  in  the  x,  y,  z  space,  we  can  consider  that  dW 
effectively  lies  in  the  viewing  plane  and  gives  the  following 
decomposition 


dH  =  dH  +  dH 

r  u 

A  A 

(along  r)  (along  u) 
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Figure  3.  Viewing  Plane  Rotation  Decomposition 


Thus  dH  determines 
u 


A 

the  direction  of  u  and 


_  h 


dR  =  r  X  ^ 


(motion^of  r 
about  u) 


u  X 


dO 


(motion^of  u 
about  r) 


(5) 


since 


r  X  dfi 

r 


A  _ 

u  X  dQ 


u 


0  . 


A 

For  example  with  u  fixed  in  inertial  space. 


du  =  u  X  dH  =  u  X  dH  =  0 

r 


so  that  u  is  parallel  to  dH,  (dH  =  .  Similarly  for  r  fixed  in 

inertial  space  ,  r  is  parallel  to  dCl  (dfi  =  dfl^)  .  In  both  cases  dH 
has  no  component  normal  to  the  viewing  plane. 

With  5(t),  the  angular  position  of  u  about  the  line  of  sight 
direction  r,  then  d?  gives  |d0^|  . 

If  the  conceptual  model  just  outlined  is  viewed  in  terms  of  a 
possible  approximation  method  for  deriving  the  motion,  not  requiring 

the  assumption  of  any  motion  model,  then  to  give  d^^  =  dH  - 

A  ^  ^ 

and  so  the  direction  of  u,  we  need  to  measure  or  approximate  r,  dw, 

and  dH  such  that 
r 
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where 


1) 


A 

r 


ft 

r 


ft  ft 

(®>  0  j  ’n)  di‘l  onrr  i  A‘d 

h 


A 

6  =  radar  elevation 

A 

0  =  radar  azimuth 

JU^  =  earth  rotation  rate 
E 

\  =  site  longitude 
T\  =  site  latitude 


2)  dfJ^  =  1  d?  I  r 


3)  dO  =  Idfi]  Iq  ;  i^  =  unit  vector  in  ^  direction 


with 

|dfi|  =  I  cos  ^  “I  lira  <  n(t  +  At)  ,  n(t)>  1  | 
^  At-0 


A 


i 


Q  “ 


lim  I 

At-0 


fo  A 

n(t  +  At)  X  n(t) 

TA -  - 

I  n(t  +  At)  X  n(t) 


j 


It  is  to  be  noted  that  the  restriction  to  a  monostatic  radar 
is  most  severe  since  at  best,  with  no  approximation,  it  requires 
the  assumption  of  a  motion  model. 

Returning  to  the  formulation  via  the  scattering  matrix,  the 
restriction  to  a  motion  model  can  be  reasonably  generalized  to 
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precession.  We  picture  this  in  Figure  4  in  terms  of  the  inertial 
frame  and  an  angular  momentum  axis  system. 

Taking  account  of  the  inertial  reference  and  angular  momentum 
systems,  Equation  (3)  can  be  written  as, 

3  3 

^  ^  CyCt)  u^(t)  -  0  (6) 

j=l  i=l 

where  Y(t)  =  tan  §(t)  is  measurable, 


u.(t)  =  u . (0 ,  0,  0^,  t) ;  j  =  1,  2,  3,  are  unknown  coordinates  of 

A  ^  ^  ^ 

u  in  the  (x ,  y,  z)  system  and  0^  the 
initial  value  of  0, 


a  .  . 
ij 


A  A 

"Hj  0(t),  0(t));  i,  j  =  1,  2,  3  are  computable 

1  j  t 

coefficients  of  the  transformation 
from  radar  to  inertial  reference 
coordinates , 


b.  . 
i-J 


b.  .0,  a); 


i,  j  =  1,  2,  3,  are  unknown  coefficients  of 
the  transformation  from  angular  momentum  to 
inertial  reference  coordinates. 


Formally  at  least  Equation  (6)  holds  for  motion  with  torque 
in  which  3  =  3(t)  ,  Cl  =  a(t)  ,  9  =  9(t)  ,  0  =  0(t)  and  if  spin  is 
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Figure  4.  General  Precessional  Treatment 


observable  Y  =  Y(t)  where  in  each  case  general  time  functions  are 

applicable.  With  torque  free  motion  3  =*  constant  and  a  »  constant 

and  6(t),  0(t)  ,  Y(t)  can  either  be  general  (e.g.  the  case  of  an  asymmetric 

top)  or,  for  the  further  specification  to  precession,  0,0,  and  ¥  are 

A 

constants.  The  formulation  of  u(t)  relative  to  3,  a,  0,  0,  0^  and  ¥ 

A 

is  straightforward  in  this  case.  That  is,  u  depends  on  six  fixed 
parameters  and  five  when,  as  here,  spin  is  assumed  unobservable.  When 
general  time  functions  are  assumed,  expansion  representations  can 
increase  the  number  of  unknown  parameters  to  an  unmanageable  extent  . 

Thus,  Equation  (6)  can  be  viewed  as  a  system,  with  time,  of  space 
curves  in  5-dimensional  space  of  the  unknown  precession  parameters 
(3,  OL,  6,  0  ,0).  Then  the  equation  is  applied  at  t  *  t,  ;  k  ^  5  for 

O  K 

solution  of  the  parameters.  This  system  is  non-linear  thus  admitting 

a  family  of  solutions  with  restrictions  imposed  by  the  allowable  range 

of  the  parameter  values. 

Generally  this  is  a  most  impractical  approach  to  solution.  Some 
help  may  be  gained  by  measuring  at  times  when  y(t)  =•  0  in  order  to 
simplify.  In  addition  a  formulation  in  terms  of  Cayley-Klein  parameters 
and  complex  direction  numbers  can  provide  some  compactness  at  least 
formally.  We  illustrate  this  by  giving  Equation  (2)  in  detail  in 
both  conventional  and  Cayley-Klein  form. 
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3  '3 


1.  ^2’  ^3’  ''®’  *^o’ 


\(t)  =  < 


i=l-  1=1 


3  3 


L  ^2’  ^3’ 


0  ,  0,  t) 


j  =  l  i=*l 


(7a) 


?(c)  =  ARG 


(13  |3  +  V  V  )  (cot  6/2  +  (-|3  V  +  V  u  ) 

ab  ab  abab 

*  *  i<  i<  i  (2)  t  ‘ ' 

(-V  |3,  +  \l  V,  )  (cot  0/2  e  )  +  (V  V  +  n  n  ) 

ab  ab  ab  ab 


(7b) 

where  cot  0/2  e  is  the  complex  direction  of  u,  initial  ^ 

the  Cayley-Klein  parameters  are 


1-1  =  cos  —  e 

a  2 


(-i/2)  (02  +  33)  ;  3^  (i/2)  (02  -  33) 

V  =  sin  -z—e 
a  2 


g  (-i/2)  (a  +  0^) 

=  cos  j  e 


Vb  •  sm  j 


j  (i/2)  (a  -  d,^) 


with 


^1’  ^2’  ^3  ”  angles  relating  inertial  system  to  radar 

system , 

0,  a,  00  =  Euler  angles  relating  inertial  system  to  momentum 

system , 


and  (  )  the  conjugate  of  (  )  . 


i0t 


The  coefficients  of  cos  0/2  e  (0  and  0t  are  the  spherical 
angles  of  the  direction)  in  the  linear  fractional  transformations 
of  (7b)  represent  the  combined  effect  of  transforming  from  the  radar 

to  the  inertial  system  via 
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a 


and  from  the  momentum  to  inertial  system  via 


by  matrix  multiplication.  The  parallelism  between  (,7a)  and  (7b)  is 
complete  if,  as  is  possible,  we  take  in  (7a) 


and 


b.  .  =  b.  .(0 ,  a,  0  ) 
ij  ij  ’  ’  o 


Uj  =  Uj  (9,  0,  0,  t) 


An  alternate  phrasing  of  (7a)  sometimes  considered  can  be 
written  as  ^ 

3 

Y(t)  =  (  -ij - 

ICijCPi’.  02' > 
j=l 


(8) 


where 


kj  .  ki  ‘  1  '  2  '  3'  '  Ij 

i=l  i=l 


k  =  1,  2 
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so  that 


P  '  =  3  '  (3,  a,  t)  ;  s  =  1,  2,  3 
s  s 


Also  5 


3  3 

cos  A(t)  =  )  )  a  .(3,  CX)TI 

Lmi  ImJ  J  1  1  J 


A  A 

where  A(t)  is  the  aspect  angle  between  u  and  r.  However,  this  relating 

A 

of  the  momentum  to  radar  systems  directly  hence  of  u  to  the  radar  in 
terms  of  aspect  angle  A(t)  is  of  little  consequence  since  the 
parameters  3,  CX  must  eventually  be  exposed  for  a  proper  solution. 

More  importantly,  all  approaches  which  attempt  to  solve  for 

the  precession  parameters  via  this  system  of  non-linear  equations 

_  A 

are  at  best  arduous.  Rather  if  this  equation,  <n,  u>  =  0  is  considered 
as  the  equation  of  a  plane  (viewing)  in  combination  with  a  priorj. 
knowledge  of  precession  mechanics  then  an  alternate  and  greatly 
simplified  approach  emerges. 

However,  before  passing  on  to  the  details  of  the  method,  a 

few  remarks  concerning  particular  cases  are  now  appropriate. 

A 

Case  (1)  u  fixed  in  space  (stabilized  inertially) . 

Then  the  intersection  of  two  viewing  planes  at  separate  times 

A 

provides  u  as  illustrated  in  Figure  5. 

Also  with  very  slow  motion  in  inertial  space,  with  no  motion 
model  assumed,  the  above  method  can  serve  as  an  approximate  solution. 
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Figure  5.  Intersecting  Viewing  Planes 


A  A 

With  u  determined  the  angular  momentum  unit  vector  direction  m  is 

T 

then  given  by  the  eigenvector  of  U  U  for  the  minimum  eigenvalue  where 


u  .  J  -  1,  2.  3;  i  -  1,  2,...,  N 


and  u^j  are  the  components  of  u(t^)  in  inertial  coordinates. 

Case  (2)  If  the  local  motion,  is  very  rapid  compared  to  other 

tiUU 

motions,  then  with  t^  -  t^  ^  =  T,  where  T  is  the  precession  period, 

A 

u(t.)  is  effectively  stationary  and  is  given  by  the  eigenvector  of  K  k 

A 

for  the  minimum  eigenvalue  where,  using  components  of  n. 


K  =  j  =  1,  2,  3;  i  =  1,  2,. 


N  +  1 


A  test  on  the  determination  of  T  can  be  made  in  terms  of  the 

A  A 

singularity  of  H,  After  u  is  obtained  then  m  can  be  obtained  from 

[u(t)  +  u(t  +  T/2) ] ,  With  all  other  precession  parameters  obtained 

A 

a  recalculation  of  T,  using  the  formulation  of  u(t)  in  terms  of 

6,0,0,  can  be  made  to  ad  lust  T  and  repeat  the  process  for  improved 
o 

estimates . 
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SECTION  III 


THE  METHOD 


Briefly  the  method  relies  on  the  correct  solution  having  equal 

values  of  actual  or  scaled  azimuth  spacing  on  the  precession  cone  of 

the  body  axis  orientation  determined  at  successive  pulse  times  that 

is  on  the  linear  relation  between  azimuth  and  time.  The  specification 

of  body  axis  orientation  is  made  as  the  intersection  of  the  half  cone 

of  precession  and  a  plane,  the  viewing  plane,  containing  the  line 

of  sight  and  the  body  axis  (the  viewing  plane  also  contains  the 

projection  of  the  body  axis  onto  the  measuring  plane) . 

More  particularly,  an  initial  trial  cone  of  precession  is  chosen 

by  specifying  with  respect  to  the  chosen  inertial  system  the  azimuth 

(cl  )  and  polar  angle  (3^)  of  the  cone  axis  as  well  as  its  half  angle 

of  precession  (0  ).  With  a  ,  3  ,  9  given,  the  intersection  of  the 
^  c  c  c  c 

cone  with  a  sequence  of  say  N  viewing  planes  provides  at  most  a 

2N 


sequence  of  2N  body  axis  directions  denoted  by  {u^}  since  each 
plane  and  the  cone  intersect  along  two  direction  lines  . 

Using  unperturbed  data  measured  at  equal  time  increments  or 


N 

appropriately  scaled  one  of  the  consistent  subsequences  of  {u.} 

^  i  =  l 

will  result  in  equal  azimuth  spacing  about  the  cone  axis  if  the  trial 

cone  parameters  ^  correct.  Equivalently  stated,  such 

a  sequence  of  azimuth  differences  will  have  zero  variance.  The 

parameters  a  ,  3  ,  and  9  on  which  the  azimuth  values  {0.}  and  so 
c  c  c  1 
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the  variance  depend  are  adjusted  to  the  optimum  values  for  minimum 
variance  which  is  zero  in  the  unperturbed  case.  Various  functional 
minimization  procedures  are  possible.  It  is  possible  to  begin  with 
a  form  of  gradient  method  known  as  direct  search,  then  as  required 
by  a  reduced  convergence  rate  to  switch  to  another  calculation,  say 
Newton-Raphson ,  near  the  vicinity  of  the  minimum. 

The  testing  of  the  method  can  be  done  with  the  aid  of  a 
radar  motion  simulation  program  developed  in  mid- 1964.  Tliis 
simulator  characterizes  motion  in  terms  of  the  earth,  orbit  and  local 
(about  object  c.m.)  motion.  The  latter  allows  for  general  torque 
free  motion  or  various  stabilized  motions.  An  acquisition  capability 
for  any  number  of  sites  and  objects  is  also  included.  It  provides 
a  scattering  matrix  radar  model  in  linear  or  circular  polarization 

i 

with  conversion,  and  with  data  either  unperturbed  or  perturbed  by 
Faraday  rotation,  signal  level  and  phase  variation  with  range  and 
various  noise  in  amplitude  and  phase.  In  addition  a  model  is  included 
for  handling  the  bistatic  case.  In  the  simulator  the  angle  §(t)  is 
obtained  as  an  Euler  angle  relating  a  body  axis  system  to  the  radar 
system  of  coordinates. 
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SECTION  IV 


FORMULATION 

From  scattering  matrix  data  the  i^^  azimuth,  of  the  body 

axis  projection  on  the  plane  normal  to  the  line  of  sight  allows  for 

a  determination  of  the  orientation  of  the  i^^  viewing  plane  in  terms 

h 


of  its  unit  normal  n,  as 

I 


AAA 
n.  -  sin  I,  -  cos  5^ 


where  x  and  y  are  unit  coordinate  vectors  in  the  radar  system, 

R  R 

A 

Written  in  the  chosen  inertial  frame,  the  components  of  n^  are 

denoted  by  n  ,  n  n.^.  The  basic  geometry  is  reviewed  in 
i  1  i2  i  J 

Figure  6, 

In  order  to  maintain  a  consistent  choice  of  test  half  cone, 

A  A  A 

<u ,  is  taken,  say,  always  positive  where  m  is  the  unit  vector  of 

the  trial  cone  axis  direction  (along  the  angular  momentum  line).  In 

addition,  it  is  useful  in  choosing  consistent  precession  angle  values, 

A  A 

0,  and  in  setting  up  tests  to  deal  with  noisy  data  to  take<n,  uf>  >0 

that  is 
3 


\  n.  .m,  =  n.-  sin  P  cos  CX  -f  n.^ 
^  LJ  J  ll  c  c  l2 

j=l 


sin  0  sin  O.  +  n..  cf.s  0 
c  c  l3  c 


=  n. „  >  0 
li 
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where 


n .  . 


A 

j  =  1,  2,  3,  are  the  coordinates  of  m  in  the  inertial  system 


A 

;  j  =  1,  2,  3,  are  the  coordinates  of  n^  in  the  test  cone  axis 

AAA 

system  (x,  y,  z) 


A  A 

z  =  m  is  related  to  the  inertial  system  by  the  polar  and  azimuth 
angles  and  OL^  respectively,  and 


0  ^  P  <  IT/2:  0  ^  a  <  2TT  . 
c  c 


If  3 

n.  .m .  <  0  , 

^  1 
j  =  l 

the  signs  of  j  =  2,  3,  are  reversed. 

In  addition  to  a  and  P  the  trial  half  cone  is  then  fully 
c  c 

specified  by  assignment  of  0^,  the  half  angle  of  precession  with 

0  <  0  ^  TT/2. 

c 

With  noisy  data,  the  intersection  of  all  viewing  planes  with 
the  correct  cone  is  not  assured.  To  allow  intersection  of  the  test 
cone  with  most  of  the  measured  viewing  planes,  if  possible,  0^  is 
chosen  so  as  to  exclude  only  a  small  percentage,  say  l-i  of  chese 
planes.  We  may  phrase  this  situation  as  follows;  let 

f2(0c)  =  {i;  ^  0^} 
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having  measure 


Further,  certain  viewing  planes  which  are  preserved  by  the 
above  selection  process  may  be  near  tangency .  Of  the  viewing  planes 
chosen  we  select  the  set  of  i  values  denoted  by  n(p,  0^)  where 


n(p,  e^)  c  n(0^) 


such  that, 

^(P,  0^)  *=  ii;  s  p  sin  0^] 


Thus,  |dg  and  p  are  chosen  to  be  less  than  but  near  1;  for 
c 

example  pig  =  ,99  and  p  =  .95, 
c 

In  order  to  visualize  the  analytic  formulation  more  clearly, 
we  show  the  geometry  in  Figure  7. 

As  noted,  for  unequal  interpulse  times,  the  values  are 

simply  scaled  for  equal  values  using  for  each  i 


K  1 
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Figure  7.  Geometry  of  the  Motion  Method 


Using  then  the  two  equations  of  the  viewing  plane  and  half 


cone  in  test  cone  coordinates  we  obtain 


A  rv/  n 

<n.  X>  =  n.,  (sin  6  cos  0)  +  n.^  (sin  9  sin  0)  +  n.„  cos  9*0 
i’  il  ^  ^  i3 


(9) 


e  =  e 


(10) 


With  9  <  tt/2,  this  establishes  (with  'I]  are  already  chosen)  consistently 

A 

the  test  half  cone.  The  sense  of  Ti  is  set  by  constraints  on  a  and  P  . 

c  c 

th  ^ 

Denoting  the  azimuth  of  the  i  determination  of  u  as  0^  in 

rs./ 

the  axis  system  (x,  y,  z)  attached  to  the  test  cone  as  previously 


described,  we  obtain  as  a  function  of  OL  ,  3  ,  and  9  only 

c  c  c 


0.  =  tan 
1 


1 

<\2] 

+  cos  ^ 

i  1 

-^3 

2  2  k  1 

(n. ,  +  n. „  )^  tan  9  / 

'  il  i2  c  ' 

(11) 


=  f.(a 

I  c 


P,.)  +  Si  (%: 


where 


n .  . 
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The  first  term  f.  measures  the  azimuth  of  the  vector  n.  while* 

1  I 

fi 

g.  measures  the  azimuth  of  u.  relative  to  n.  with  both  measu renen ts 

L  1  L 

in  the  x,  y  plane  of  the  trial  cone  frame. 

Since,  with  n^^  ^  ^  » 


-n 


i3 


(n. -  +  n. ^  )  tan  ^ 

1 1  i2  c 


then , 


;/2  ^  g.  ^  3r/2 


which  allows  two  solutions  for  u^  symmetric  about  g^  -  We 

denote  these  solutions  as  g.,  and  where  the  consistent  selection 

il  °i2 

of  f  >  g  ^nd  so  0  values  is  based  on  sign  changes  in  n^^  ,  the  possible 
monotonicity  of  0  ,  and  consistent  (  A0)  values. 


TT/2  ^  ^  TT  ,  11  ^  g.2  iU/2  . 


Using 


0..  =  f.  +  B-k!  k  =  1,  2, 

X  K  1  X  K 


we  calculate 


j-i 
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whe  re  , 


j  =  MIN  [  i '  ;  i '  >  i] 

i' 


for  all  i  allowed  by  the  selection  rules  with  k  =  1,  2  and  take  j-i 
calculations  for  each  i,  j  pair. 

The  azimuth  relationships  are  shown  in  Figure  8. 

Thus  using  the  criterion  of  equal  A0  values,  a  function  is 
constructed  which  for  the  proper  g  set  and  correct  precession 
parameters  is  minimum. 


M 

Denoting  the  set  of  such  differences  as  [A0  ,  3  5  ^  =  1 ,  2 


mk  m=l 


then  we  calculate  for  k  =  1,  2; 


A0  the  mean 

mk 


m=l 


and 


2 


M 


a 


A0, 


M 


m=l 


—  2 

(A0  ,  -  A0  )  the  variance 

mk  k 


where 


2 


k 


A0, 


k 
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Figure  8.  Azimuth  Relations 


Notice,  the  formulation  is  derived  in  terms  of  just  three 


precession  parameters  so  that  we  seek 


*  *  Vc 

0  ,  a  ,  0 

c  c  ’  c 


to  give 


C-  C-  >  =  “c 


(12) 


A0 

If  we  take  0  '  =  0  +  A0  with 

c  c  c  y 


«1  and,  calculate 


0  ’  ^c  ’  ^c^^’  ^  ^  after  applying  the  selection  rules 

k 

based  on  Idg  and  p,  and  proceed  similarly  for 
c 

(0-  '  ,  P  ,  0  )  and  (a  ,  0  '  ,  0  ) 

A0,  ^  c  ’  c’  C  A0,  ^  c’  C  ’  c 

;k  k 


also  including  in  these  two  cases  the  proper  initial  adjustment  of 

2 

the  sign  of  then  the  function  G  can  be  altered  by  simultaneous 

changes  in  coordinates  j  9^)  along  a  direction  approximating 

steepest  descent.  Such  a  method  generating  so  called  normal  equations 

d  ^ 

is  often  inefficient  in  computing.  In  these  methods  y  evaluated 
or  approximated  for  (  )  equal  to  OL^ ,  or  0^  gives  the  new  value 
directly  as , 


Aa  +  AP  +  Ij-  ^0  +  (a  ,  0  ,  0  ) 

oa  c  op  c  00  c  A0  c  ’  c  c 
c  c  c  k 


k  =  1,  2 
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The  actual  algorithm  used  initially  will  follow  a  method  called 

direct  search.  Adaptation  and  revisions  of  this  approach  have 

been  considered  in  internal  MITRE  reports.  This  method  uses  two  basic 

approach  moves  called  exploratory  and  pattern.  Each  pattern  move  follows 

a  successful  exploratory  move  by  duplicating  this  movement  along  the 

★ 

direction  of  change  established  by  that  exploratory  move. 

If  an  exploratory  move  cannot  effect  a  decrease  in  the  function 
2 

to  be  minimized,  here  O  ,  the  coordinate  position  is  returned  to  that 
effected  by  the  previous  exploratory  move.  Then  a  new  exploration 
is  made  which  if  unsuccessful  results  in  a  reduction  of  coordinate 
step  size  and  another  exploration  etc.  In  each  exploration  the 
coordinate  size  is  fixed  and  its  sign  is  chosen  to  cause  a  function 
decrease  if  possible.  Each  coordinate  is  treated  sequentially 
starting  with  the  function  value  established  by  the  previous 
coordinate  change . 

Each  pattern  move  whether  successful  or  not  is  followed  by  an 
exploratory  move  from  a  new  base  point  upon  success  or  the  former 
base  point  on  failure. 

It  appears  that  other  useful  additions  to  this  process  might 
check  on  the  rate  of  change  with  respect  to  the  coordinates,  in  the  manner 

"^This  is  a  variation  of  the  usual  approach  of  duplicating  both  the 
previous  pattern  and  exploratory  moves.  For  a  sequence  of  successful 
pattern  moves  this  could  amount  to  duplicating  the  effect  of  all 
previous  exploratory  moves.  A  preferred  approach  would  seem  to  be 
one  which  weights  recent  exploratory  moves  more  favorably. 
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outlined  above  under  the  discussion  of  partial  derivatives  to 
establish  initial  coordinate  step  size  ratios  as  well  as  to  determine 
when  to  switch  to  another  method  if  the  step  sizes  are  too  small  in 
the  vicinity  of  solution. 

In  summary  the  motion  method  is  as  follows: 

2 

1.  Select  initial  3  ,  a  ,  0  values  and  compute  (3  >  ^  >  9  ) 

C  C  C  tC  c  c  c 

applying  selection  rules.  Other  calculations,  such  as  intersecting 

planes,  provide  aid  in  choosing  initial  values. 

*  "k  Vc 

2.  Apply  minimization  procedure  to  give  correct  ° ^  ^ 
values  for  k  =  1,  2,  using  for  example  modified  direct  search  and 
Newton-Raphson  methods.  Redundancy  in  the  data  with  the  correct  k 
value  can  be  used  to  improve  estimates, 

3.  Establish  0,  0^  directly  from  0^  values  and  the  optimum 

^  ^  -k 

set  3  ,  CX  ,  0 

c  c  c 

Obvious  extensions  to  the  bistatic  case,  using  the  simulator 
again  as  a  test  mechanism,  can  be  made.  The  testing  procedure  will 
emphasize  unperturbed  data  to  test  the  accuracy  and  usefulness  of 
the  basic  method.  This  will  be  followed  by  tests  with  **noisy  data’’ 
to  evaluate  the  capability  for  practical  operation. 
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